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Abstract 

Square grid circle patterns with prescribed intersection angles, mimicking holomorphic maps 
and log(z) are studied. It is shown that the corresponding circle patterns are embedded 
and described by special separatrix solutions of discrete Painleve and Riccati equations. 
General solution of this Riccati equation is expressed in terms of the hypergeometric function. 
Global properties of these solutions, as well as of the discrete and log(z), are established. 


1 Introduction 

The theory of circle patterns is a rich fascinating area having its origin in classical theory of 
circle packings. Its fast development in recent years is caused by mutual influence and interplay 
of ideas and concepts from discrete geometry, complex analysis and the theory of integrable 
systems. 

The progress in this area was initiated by Thurston’s idea mm about approximating 
the Riemann mapping by circle packings. Classical circle packings comprised of disjoint open 
disks were later generalized to circle patterns where the disks may overlap (see for example 
M)- Different underlying combinatorics were considered. Schramm introduced a class of circle 
patterns with the combinatorics of the square grid HU; hexagonal circle patterns were studied 
in and [Jj. 

The striking analogy between circle patterns and the classical analytic function theory is 
underlined by such facts as the uniformization theorem concerning circle packing realizations of 
cell complexes of a prescribed combinatorics jlj, discrete maximum principle, Schwarz’s lemma 
uni and rigidity properties discrete Dirichlet principle EH. 

The convergence of discrete conformal maps represented by circle packings was proven by 
Rodin and Sullivan m- For a prescribed regular combinatorics this result was refined. He and 
Schramm m showed that for hexagonal packings the convergence is C°°. The uniform conver¬ 
gence for circle patterns with the combinatorics of the square grid and orthogonal neighboring 
circles was established by Schramm EH- 

Approximation issue naturally leads to the question about analogs to standard holomorphic 
functions. Computer experiments give evidence for their existence mm however not very 
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much is known. For circle packings with the hexagonal combinatorics the only explicitly de¬ 
scribed examples are Doyle spirals mM which are discrete analogs of exponential maps and 
conformally symmetric packings, which are analogs of a quotient of Airy functions [2]. For 
patterns with overlapping circles more explicit examples are known: discrete versions of exp(z), 
erf(z) |21 j . z"^, log( 2 :) |2j are constructed for patterns with underlying combinatorics of the square 
grid; z'^, log( 2 ;) are also described for hexagonal patterns [HI, [li¬ 
lt turned out that an effective approach to the description of circle patterns with overlapping 
circles is given by the theory of integrable systems (see 0,0,0). For example, Schramm’s circle 
patterns are governed by a difference equation which is the stationary Hirota equation m (see 
m for a survey). This approach proved to be especially useful for the construction of discrete 
z'^ and log( 2 ;) in 0,0,0 with the aid of some isomonodromy problem. Another connection 
with the theory of discrete integrable equations was revealed in 0,0: embedded circle patterns 
are described by special solutions of discrete Painleve II equations, thus giving geometrical 
interpretation thereof. 

This research was motivated by the attempt to carry the results of 0,0 over square grid 
circle patterns with prescribed intersection angles giving Schramm patterns as a special case. 
Namely, we prove that such circle patterns mimicking z"^ and log(z) are embedded. This turned 
out to be not straightforward and lead to asymptotical analysis of solutions to discrete Riceati 
equation. As the Riccati differential equation is known as the only equation of the first order 
possessing Panleve property we are tempted to conclude that circle patterns are described by 
discrete equations with Painleve property though there is no satisfactory generalization thereof 
on discrete equations. 

We use the following definition for square grid circle patterns, which is slightly modified version 
of one from m- 

Definition 1 Let G be a subgraph of the 1-skeleton of the cell complex with vertiees Z + iZ = Z^ 
and 0 < a < TT. Square grid circle pattern for G with intersection angles a is an indexed 
eollection of circles on the eomplex plane 

{C,:ze V{G), G, e C} 


that satisfy: 

1) if z, z + i G V{G) then the intersection angle of Gz and Gz+i is a , 

Sj if z, z -h 1 G V(G) then the intersection angle of Cz and Gz+i is tt — a, 

3) if z,z-\-l-\-iGV (G) (or z, z — 1 -\- i G V {G)) then the disks, defined by Cz and Cz+i+i (Cz 
and Cz-i+i respectively) are tangent and disjoint, 

4) if z,zi,Z 2 G V{G), \zi — Z 2 \ = \/ 2 , jz — zil = jz — Z 2 l = 1 (i.e. Czi,Cz 2 are tangent and Cz 
intersects Czi and and Z 2 = z + i{zi — z) (i.e. Z 2 is one step counterclockwise from zi), 
then the circular order of the triplet of points Cz H Cz^ — Cz 2 ,Czi n Cz 2 ,Cz H Cz 2 — Cz^ agrees 
with the orientation of Cz. 

The intersection angle is the angle at the corner of disc intersection domain (Fig. 0). 

To visualize the analogy between Schramm’s circle patterns and conformal maps, consider regular 
patterns composed of unit circles and suppose that the radii are being deformed to preserve the 
intersection angles of neighboring circles and the tangency of half-neighboring ones. Discrete 
maps taking the intersection points and the centers of the unit circles of the standard regular 
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Figure 1 : Circles intersection angle. 



Figure 2 : Schramm type circle patterns as a discrete conformal map. The discrete version of 
the holomorphic mapping . The case tana = 3 . 


patterns to the respective points of the deformed patterns mimic classical holomorphic functions, 
the deformed radii being analog of \f'{z)\ (see Fig. [^. 

In section | 2 l we give dehnition of discrete Z'^ as a solution to some integrable equation 
subjected to a non-autonomous constraint. Its geometrical properties of immersion and embed¬ 
dedness are expressed in terms of solution for radii of corresponding circle patterns. 

For this solution to be positive it is necessary that some discrete Riccati equation has positive 
solution. This equation is studied in section |31 where its general solution is expressed via the 
hypergeometric function. 

Section | 1 ] completes the proof of embeddednes, discrete equations of Painleve type being the 
main tool. Possible generalizations for non-regular combinatorics and for non-circular patterns 
are discussed in section |H1 


2 Discrete and square grid circle patterns of Schramm type. 

Definition 2 Discrete map Z'^, 0 < 7 < 2 is the solution f : C of 

q{fn ,m;/n+l,m;/n+l,m+l;/rijm+l) — ^ ( 1 ) 

n r\ - fn ,m) {fn.m - fn-l - fn ,m) {fn.m - fn ,m— 1 ) /rx\ 

lfn,m = 2n - - - - -h 2m - - - - - , (2) 


/n+l,m fn—l,r 


- fn ,m—1 
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with 0 < a < vr and the initial conditions 


/i,o 


1 , /o,i = e 


— 0*7“ 


(3) 


where q stands for cross-ratio of elementary quadrilaterals: 


9(/i,/2,/3,/4) 


(/l-/2)(/3-/4) 

(/2-/3)(/4-/l)’ 


Constraint 0 was obtained from some isomonodromy problem of Lax representation of o, 
found for 7 = 1 in na (see also 1210) The definition of can be justified by the following 
properties: 

• if one thinks of / as defined on the vertices of the cell complex with diamond-shaped faces 
(see FigEI) then o means that / respects the cross-ratios of the faces and therefore is 
’’locally conformal”, 

• asymptotics of © as n,m —> 00 suggests that / approximates z'^. 


Remark. Equation with a = tt!2 was used in [H] to define discrete conformal maps. 

The motivation was that / maps vertices of squares into vertices of the ’’conformal squares”. 
Consider the surface glued of these conformal squares along the corresponding edges. This 
surface is locally flat but can have cone-like singularities in vertices. If the map is an immersion 
then the corresponding surface does not have such singularities. Therefore it is more consistent 
to define as discrete conformal an immersion map on the vertices of cell decomposition of C 
which preserves cross-ratios of its faces. 


Proposition 1 Q© Constraint m is compatible with m- 


Compatibility is understood as a solvability of some Cauchy problem. In particular a solution 
to 0,0 in the subset is uniquely determined by its values fi^, /oq. Indeed, constraint 0 
gives /o,o = 0 and defines / along the coordinate axis (n, 0 ), ( 0 , m) as a second-order difference 
equation. Then all other fk,m with {k, m) G are calculated via cross-ratios ©. 

In this paper the following initial conditions for 0,0 are considered: 

/i.o = 1, /o,i = (4) 

with real (5. The solution / dehnes a circle pattern with square grid combinatorics. For a = 7112 
such circle patterns were introduced by Schramm in HU. For any 0 < a < vr we obtain a natural 
generalization of Schramm circle patterns. 

In what follows we say that the triangle (.zi, 22 , 2 : 3 ) has positive (negative) orientation if 


^3 - zi 
Z2 - Zl 


Z3 - Zl 
Z2 - Zl 


with 0 < (f < TT 


( —TT < (/> < 0 ). 


Lemma 1 Let q(zi, 22 , 23 , 24 ) = e 0 < a < vr. 

• // l^i — 22 I = | 2 i — 24 I and the triangle ( 24 , 22 , 24 ) has positive orientation then 123 — 22 
123 — 241 and the angle between [ 2 : 1 , 22 ] and [ 22 , 23 ] is (vr — a). 
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• If \zi — Z2\ = \zi — z^l and the triangle {zi, Z2, z^) has negative orientation then \z^ — Z2\ = 
1^3 — Z4I and the angle between [-2;i,-Z2] and [-2;2,-Z3] is a. 

• If the angle between [-21,-22] and [-21,24] is a and the triangle (21,22,2^4) has positive orien¬ 
tation then I23 — 22I = |2i — 22I and [23 — 24] = [24 — 211. 

• If the angle between [21,22] and [21,24] is (tt — a) and the triangle (21,22,24) has negative 
orientation then I23 — 22] = |2i — 22] and [23 — 24] = [24 — 2i|. 

Proof: straightforward. 

Proposition 2 All the elementary quadrilaterals (/n,m,/n+i,m,/n+i,m+i,/n,m+i) for the solu¬ 
tion of WM with initial m are of kite form: all edges at the vertex fn,m with n-\-m = 0 (mod 2) 
are of the same length. Moreover, each elementary quadrilateral has one of the forms enumerated 
in lemmaUi 

Proof: Given initial /o,i and fifi constraint (O gives fn,o and /o,m for all n, m > 1 . It is easy to 
check that / has the following equidistant property: 

f 2 n ,0 /2n—1,0 — /2n+l,0 /2n,0) fo, 2 m fo, 2 m—l — /o,2m+l fo, 2 m (5) 

for any n > 1 , m > 1 . Lemma^and |/i^o ~ /o,o| = |/o,i ~ /o,o| allows one to apply induction in 
n, m, starting with n = 0, m = 0. 

Proposition O implies that for n + m = 0 (mod 2 ) the points fn±i,m,fn,m±i he on the circle 
with the center at fn,m- For the most (3 (namely for f 3 ^ a) the behavior of thus obtained circle 
pattern is quite irregular: inner parts of different elementary quadrilaterals intersect. 

Definition 3 A discrete map fn,m is called an immersion if inner parts of adjacent elementary 
quadrilaterals {fn,m,: fn+i,m,i fn+i,m+ii fn,m+i) ore disjoint. 

Consider the sublattice {n, m : n + m = 0 (mod 2 )} and denote by V its quadrant 

Y = {z = N + iM : N,M eZ^,M > |AI|}, 

where 

N = {n — m)/ 2 , M = (n + m)f 2 . 

We use complex labels z = N + iM for this sublattice. Denote by C'(2) the circle of the radius 

Hz — \ fn,m /n+l,m| — \ fn,m /n,m+l| — \ fn,m fn—l,m\ — \ fn,m /n,m—ij (6) 

with the center at fN+M,M-N = fn,m- 

Let {C(2)}, 2 G V be a square grid circle pattern on the complex plane. Define fn,m '■ 
—> C as follows: 

a) if n + m = 0 (mod 2) then fn,m is the center of C'(^^-^ + 

b) if n + m = 1 (mod 2) then— C(M^+ nC(M±i+f^^±f±i) = G(M±i + 

i )nC( +i ). Since all elementary quadrilaterals {fn,m, /n+i,m, /n+i,m+i, /n,m+i) 
are of kite form equation m is satisfied automatically. In what follows the function fn,rm defined 
as above by a) and b) is called a discrete map corresponding to the circle pattern {C(2)} . 
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Proposition 3 Let the solution f of mM with initial 0 is an immersion, then R{z) defined 
by satisfies the following equations: 



(TV + M){R:z^i + i?2+i)(T?2 — Rz+iRz-i + cosaRz{Rz-i — Rz+i)) + 

(M — N){Rz-i + Rz+i){Rz — Rz+iRz+i + cos aRz{Rz+i — Rz+i)) = 0 , 


( 8 ) 


for z E Yrint := V\{±iV + iN\N E N}. 

Conversely let R{z) : V —> R+ satisfy 0 for z gYi and for z E Yrint- Then R{z) define a 
square grid circle patterns with intersection angles a, the corresponding discrete map fn,m is an 
immersion and satisfies GDM- 

Proof: Circle pattern is immersed if and only if all triangles {fn,m, fn+i,m, fn,m+i) of elemen¬ 
tary quadrilaterals of the map fn,m have the same orientation (for brevity we call it orien¬ 
tation of quadrilaterals). Suppose that the quadrilateral {fo,o, fi,o, fi,i, fo,i) has positive ori¬ 
entation. Let the circle pattern fn,m be an immersion. For n -|- m = 1 (mod 2 ) points 
fn,m, fn-i,m+i-, fn-2,m, fn-i,m-i he on circle with the center at fn-i,m and radius Rz, where 
z = {n — m — l )/2 + i{n -|- m — l )/2 (See the left part of Fig. O}. Using o one can compute 
fn,m+i and fn,m-i- Lemma [T] and proposition 12 imply that fn+i,m is in line with fn-i,mjn,m 
and that the points fn,m+i-,fn,m-,fn,m-i are also collinear. Denote by Rz+i-,Rz+i the radii of 
the circle at /„ m-i and /„ m+i respectively. Define Rz+i+t = Rz . If Q is satisfied 

at (n — 1, m) then 0 at (n, m) is equivalent to 0, Rz+i+i being positive iff the quadrilaterals 
{fn,m: fn+l,m: fn+l,m+l: fn,m+l') and {fn,m—l:fn+l,m—l:fn+l,m:fn,m) have positive orientation. 

Similarly starting with 0 at {n,m — 1), where n -|- m = 0 (mod 2) (see the right part of 
Figini one can determine evolution of the cross-like figure formed by fn,m-i, fn+i,m-i, fn,m, 
/n—l,m—1) fn,m—2 iuto fn+l,mi fn+2,mi fn+l,m+l} fn,mt fn+l,m,— l- Equation (|2) at (tI-|- l,?Ti) is 
equivalent to 0 and Q at z = (n — m)/2 -|- i{n + m)/2. Rz+i is positive only for immersed 
circle pattern. 

Now let Rz be some positive solution to 0,®. We can rescale it so that Rq = 1 . This 
solution is completely defined by Ro,Ri. Consider solution fn,m of 0,0 with initial data 
0 where f 3 is chosen so that the quadrilateral (/o,05/i,0)/i,i)/o,i) has positive orientation and 
satisfies the conditions Rq = 1 = |/o,o — /i,o| and Ri = |/i^i — /i,o|- The map fn,m defines circle 
pattern due to proposition |2 It can be uniquely computed from 0,®. To this end one have to 
resolve 0 with respect to Rz+i+i and use it to find fn+i m from Rz+i+i = ) and 

to resolve ® for Rz+j to find fn+i,m+i from Rz+i = Rz+i • One can reverse 

the argument used in derivation of 0 ,® to show that / satisfies 0 , 0 . Moreover, since Rz 
is positive, at each step we get positively orientated quadrilaterals. Q.E.D. 

Note that initial data @ for fn,m imply initial data for R^: 



(9) 


Definition 4 A discrete map fn,m is called embedded if inner parts of different elementary 
quadrilaterals {fn,m, fn+i,m, fn+i,m+i, fn,m+i) do not intersect. 
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Figure 3 : Kite-quadrilaterals of circle pattern. 


Proposition 4 If for a solution Rz of with 7/1 and initial conditions m holds true: 

Rz >0, (7 - l){Rl - Rz+iRz-i + cos aRz{Rz-i - Rz+i)) > 0 (10) 

in 'Vint, then the corresponding discrete map is embedded. 

Proof: Since R{z) > 0 the corresponding discrete map is an immersion. Consider piecewise 
linear curve r„ formed by segments [fn,m, fn,m+i], where n>0, 0<m<n — 1 and the vector 
'Vn{m) = {fn,mfn,m+i) along this curve. Due to Proposition lathis vector rotates only in vertices 
with re -|- m = 0 (mod 2 ) as rre increases along the curve. The sign of the rotation angle 9 n{m), 
where —vr < 0 n{m) < vr, 0 < rre < re is defined by the sign of expression 

Rz Rz+lRz—i T cos CxRz{Rz—i Rz+l) 

(note that there is no rotation if this expression vanishes), where z = {n — rre)/2 -|- i{n + rre)/2 
is a label for the circle with the center in fn,m- If re -|- rre = 1 (mod 2 ) define 0 n{m) = 0 . 

Now the theorem hypothesis and equation (jS)) imply that the vector Vn{m) rotates with 
increasing rre in the same direction for all re, and namely, clockwise for 7 < 1 and counterclockwise 
for 7 > 1. Due to the following Lemma it is sufficient to prove m only for 1 < 7 < 2. 

Lemma 2 If Rz is a solution of 0,® for 7 then l/Rz is a solution of for 7 = 2 — 7. 

Consider the sector B := {z = rel^ : r > 0 , 0 < (/? < 7a/2}. The terminal points of the 
curves T^ he on the sector border. 

Lemma 3 The eurve r„ has no self-interseetion and lies in the seetor B. 

Lemma is proved by induction. For re = 1 it is obviously true since the curve Fi is a line segment. 
Let it be true for re > 1 . Note that for immersed / the first and the last segments of F„_|_i lie in 
B therefore the vector v„_|_i(rre) rotates counterclockwise around convex domain of B confined 
by F„ and can not make a full circle. Lemma is proved. 
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Each curve E^ cuts the sector B into a finite part and an infinite part. Since the curve r„ 
is convex and the borders of all elementary quadrilaterals {fn,m, fn+i,m, fn+i,m+i, fn,m+i) for 
imbedded Z'^ have the positive orientation the segments of the curve r„_|_i lie in the infinite 
part. Now the induction in n completes the proof of Proposition | 1 ] for 1 < 7 < 2 . The proof for 
0 < 7 < 1 is similar. Q.E.D. 

One can compute q from m and obtain the asymptotics q ~ c{'^)irB from Stirling 
formula for large x: 

r(x) ~ . (11) 

Numerical experiments support the following conjecture. 

Conjecture. 


3 Discrete Riccati equation and hyper geometric functions 

Let Tn and Rn be radii of the circles with the centers at f2n,o, /2n+i,i respectively (see Fig ID). 



Constraint © and property © gives 

Tn+l = 


2n + 7 


2(n + l)-7 

From elementary geometric considerations one gets 

^n+l Rn COS Oi 
■^n+1 — ^n+1 


( 12 ) 


Define 


Rn ^n+1 COS CX 
Rn , . 2n + 7 


and denote t = cos a for brevity. Now the equation for radii i?, r takes the form: 


Onil) - tpn 
Pn -tgnilY 


( 13 ) 


Remark. Equation ( 1131 ) is a discrete version of Riccati equation. This title is motivated by 
the following properties: 








• cross-ratio of each four-tuple of its solutions is constant as Pn+i is Mobius transform of Pn, 

• general solution is expressed in terms of solution of some linear equation (see below this 
linearisation). 

Below we call m d-Riccati equation. 


Proposition 5 Solution of discrete Riccati equation is positive for n >0 iff 


sin 


'ya 


Po = 


sm 


( 2 - 7 )q 


(14) 


Proof is based on the closed form of the general solntion of d-Riccati linearisation. It is linearised 
by the standard Ansatz 

(15) 


Pn = +tgn{'y) 


Vn 


which transforms it into 


Vn+2 + t{gn+l{l) + l)yn+l + {V - l)gn{l)yn = 0. (16) 

One can guess that there is only one initial value pq giving positive d-Riccati solution from the 
following consideration: gni'j) —> 1 as n —> oo, and the general solution of (II 61 ) with limit values 
of coefficients is = ci(— 1)”(1 -|- t)” -|- C2(l — t)"'. So pn = + ignij) —1 for ci / 0. 

However ci,C2 defines only asymptotics of a solution. To relate it to initial values one needs 
some kind of connection formulas. Fortunately it is possible to find the general solntion to (tTHl) . 


Proposition 6 The general solution to USD is 


yn = (ciAr--v^F(4^. 1^, 1 + 




2 ’ 2 ’ 2 


(17) 


„+ i _.^/2 3-7 7 - 1 1 

+C2X2 Py n :^^,7 ;-n,Z 2 ) 


2 ’ 2 ’ 2 

where Xi = —t— 1 , A2 = 1 —t, zi = (t— 1 )/ 2 , Z2 = —(l-|-t )/2 and F stands for the hypergeometric 
function. 

Proof: Solution was found by slightly modified symbolic method (see cni for method description). 
Substitution 

yn = X = n + l- 7/2 ( 18 ) 

transforms m into 

A^(x -h l)xUa;+2 + 2 t(x -h ^ ^ )xUx+l + {V - l){x + - l){x l)Ux = 0 . ( 19 ) 

We are looking for solution in the form 


OO 

'^x — ^ ^ ^m'^x,7n 

m=—oo 


(20) 
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where Vx^m satisfies 




( 21 ) 


Remark. Note that the label m in CT is running by step 1 but is not necessary integer 
therefore Vx^m is a straightforward generalization of = {x + m — l)(x + m — 2 )...{x + l)x 
playing the role of x™' in the calculus of finite differences. General solution to m is expressed 
in terms of F-function: 


r(x + m) 

Vx,m — C ^77 7 

r(x) 


( 22 ) 


Stirling formula (tTTl) for large x gives the asymptotics for Vx^m 


Vx^m — ex™ for X —> 00. 


( 23 ) 


Substituting (HB into P|l. making use of m and collecting similar terms one gets the 
following equation for coefficients: 


(A^ + 2 tA + t^-l)am.-2 + 2 ( ^ -m)(tA + f^-l)am-i + (t^-l)(l-m)(7-l-m)am = 0 . ( 24 ) 


Choice Ai = —t — 1 or A2 = 1 — t kills the term with 0^-2- To make series (EOl) convergent we 
can use the freedom in m to truncate on one side. The choice mGZormGy + Z leads to 
divergent series. For m G 4 ^ _|_ z equation (EH) gives a7+i_|_^ = 0 for all non-negative integer k 
and 


0‘l+l_k-l 


1 -t^ 
tX +— 1 


(fc_I^)(fc_l + 7 ^) 

2 k 


0-2+1 h 


( 25 ) 


where A = Ai, A2. Substitution of solution of this recurrent relation in terms of F-functions and 
(1^ yields 




1 - 


F(fc - ^)F(fc - 1 + ^)F(x + ^ - fc) 


^^\2{tX + t‘^ -1) 


F(A:)F(x) 


( 26 ) 


Lemma 4 For both X = —t — 1,1 — t series converges for all x. 

Proof of Lemma^ Since 2; = 2{t\+p-i) ~ ~ + ^)/2 Ai, A2 respectively and 

t = cos a < 1 the convergence of (1^ depends on the behavior of —— 2 ^+^2 )r(^+ 2 — 

Stirling formula CD) ensures that this expression is bounded by for some c an efix,^) 

which gives convergence. 

Series (EH is expressed in terms of hypergeometric functions: 




7-1 


>z) 


where 


F(l- 


7 — 1 7 — 1 


,1 - x + 


7-1 


,z) = l- 


(l- 7 ^)( 7 ^) 

\l-{x + ^)) 


+ ... + 


( 27 ) 


10 

















-.k L 


(1 - ^)(2 - 


+ ... 


(l_(^7^))...(fe_(a:+7^)) 


Here we use the standard designation F(a, b, c, z) for hyper geometric function as a holomorphic 
at z = 0 solution for equation 


z(l — z)Fzz + [c — (o + 6 + l)z\Fz — abF = 0 . 


( 28 ) 


Now we can complete the proof of proposition IHl Due to linearity general solution of is 
given by superposition of any two linear independent solutions. As was shown each summond 
in m satisfies the equation To finish the proof of Proposition El one has to show that 

the particular solutions with ci = 0, C2 / 0 and ci / 0, C2 = 0 are linearly independent, which 
follows from Lemma 0 


Lemma 5 As n —> oo function m has the asymtotics 

~ (n + 1 - 7/2)^(ciAi+^“^/" + C2A2+^"^/") ( 29 ) 

Proof: For n ^ oo series representation ( 1 ^ 71 ) gives F{^^, A — n,zi) ~ 1 . Stirling formula 
(fTT|) defines asymptotics of the factor (see (|2^ ).Q.E.D. 

Proof of Proposition\^ For positive pn it is necessary that c\ = 0 : it follows from asymptotics 
(El substituted into m- Let us define 


s{z) = 1 + z 


(1 7 ^)( 7 ^) 


+ ... +P 




..(l 7 ^)( 7 ^)(fe l + 7 ^) 


fc !(/ c - 2)...2 


( 30 ) 


It is the hypergeometric function F{^-^, ^ — n,z) with n = 0. A straightforward manipu¬ 

lation with series shows that 


Po = 1 + 


2 ( 7 - 1 ) 

2-7 


z + 


4 z{z — 1 ) s'{z) 
2 — 7 s{z) 


( 31 ) 


where z = 1 ^. Note that po as a function of z satisfies some ordinary differential equation 
of first order since satisfies Riccati equation obtained by reduction of (ITHl) . Computation 

sin ^ . . 

shows that —satisfies the same ODE. Since both expression (13111 and (II 411 are equal to 1 
2 

for z = 0 they coincide everywhere. Q.E.D. 

Corollary 1 If there exists immersed fn,m satisfying it is defined by initial data m- 


4 Embedded circle patterns and discrete Painleve equations 


Let Rz be a solution of o and 0 with initial condition 0. For .2; G ^int define Pn^m — Pz — 
Qn,m = Qz = rewritten as follows 


{N - M)Qn,m{'^ + Pn,m)(.Qn,m - Pn,m cos a) - {M + N)P]sf^MSN,M 
QnmK.^ + ^)Sn,m - {M - A^)(l -I- PN,M)i^N,M - Qn,m cos a)] 


( 32 ) 
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( 33 ) 


p _ (2M + ^)Pn,m + (2iV + 7)QjV,MQjV,M+i 

N,M+i + 1) - ^)Pm,m + ( 2 (M + 1) - i)Qn,mQn,m+^ ’ 

where 

5'Af,M = QAr,M “ -P/V,M + Qn,m{'^ — Pn,m) “• 

Property cm) for (IH^ . (I 33 I) reads as 

(7 “ l)(QAr,M “ Pn,m + Qn,A di^ ~ Pn,m) ^os a) > 0 , Qn,m > 0 , Pn,m > 0 . ( 34 ) 

Equations (jSlI),(Em can be considered as a dynamical system for variable M, where due to m 


Pn+i,n 


2 N + 7 
2(Ar + i)-7- 


( 35 ) 


Proposition 7 For each > 0 there exists qjy > 0 such that the solution of system 
subjected to \d 5 \) and Qn+i,n = < 1 N has property \ 34 \) for all M > N. 

Proof: Due to Lemma |21 it is sufficient to prove (13411 only for 0 < 7 < 1 . Define real function 
F{P) on R+ implicitly hy F^ — P + F {1 — P) cos a = 0 for 0 < P < 1 and by F{P) = 1 for 

1 < P. 


Q 


D„ 


-1 -0.5 


0-5 1 1.5 


Dh 


Figure 5 : The case cos a = — 1 / 2 . 


Designate 


:= {{P, Q):P> 0 ,Q> F{P)}, := {(P, Q) : Q < 0 }, 

Do := {(P, Q) : P > 0,0 < Q < P(P)}, Df := {(P, Q) : P < 0 , Q > 0 } 

as in FigEl 

System (j^ . d^ defines the map : {Pn,m,Qn,m) {Pn,m+i-,Qn,m+^)- 

has the following properties: 

• it is a continuous map on Dq. Values of ‘h^r(M) on the border of Dq are defined by 
continuity in RP^. 

• For {Pn,MiQn,m) S Po holds true ^n{M,Pn,m,Qn,m) G Dq U P„ U Dd, i.e. (P, Q) can 
not jump in one step from Dq into Df. 


12 






Consider the solution to (ESI,® with initial conditions on the segment S{N) determined by 
(ESI) and Qn,n+i = q, where 0 < q < F{Pn+i,n) then {Pn,m^Qn,m) = 

Define Sn{N) = {q : {PN,M{q),QN,M{q)) £ Dq\/N < M < n}. Then Sn{N) is a closed set 
as ^n{M) is continuous on Dq and {P,Q) can not jump in one step from Dq on the half-line 
P = 0 ,Q > 0 . As a closed subset of a segment S{N) the set Sn{N) is a collection of disjoint 
segments {Sl^{N)}. 

Lemma 6 There exists sequence {Sn'^\N)}n>N such that: 

• Sn"‘\N) is mapped by {PN,n{q), QN,n{q)) onto some curve Cn £ Dq with one terminal point 
on the curve Q = F{P) and the other on the line Q = 0 . 

• C Sn^'"\N). 

Lemma is proved by induction. For n = N + 1 it is trivial. Let it holds true for n. maps 

{P, Q) £ Dq with Q = F{P) into Du and (P, 0) with P > 0 into D^. Thus ^n{M) maps c„ into 
some curve c^+i lying in Dq U Du U D^ with one terminal point in Du and with the other in D^. 
Therefore at least one of the connected components of c^+i H Dq has its terminal points on the 
border of Dq as stated by Lemma. 

As the segments of i^)} constructed in LemmaElare nonempty there exits q^ £ Sn{N) 
for all n> N. For this qN holds true: {PN,M{q),QN,M{q)) £ Po\ 9 Po- If {PN,M{q),QN,MiQ)) ^ 
ODq then {P]sf^M+i{q),QN,M+ii^)) would jump out of Dq. Q.E.D. 

Theorem 1 The discrete map , 0 < 7 < 2 is embedded. 


Proof: Proposition 13 ensures that for any > 0 there exist Rjsf+iN, RN+i{N+i) such that the 
solution to mM with these initial values is positive for z = N+iM, M > N. Then equation 0 
defines recursively Rk+ik > 0 , K > N and implies that Rn+im > 0 far all n, m : n > A^, m > n. 
Therefore asymptotics E 3 i implies that Rn+in/RN+ i{N+i) = Qn,n+i is exactly as defined by 
initial condition for Z'^. Proposition | 1 ] completes the proof. 

Remark 1. For A" = 0 system EUiESl for Qn,m, Pn,m reduces to the special case of discrete 
Painleve equation (the case a = nj 2 was studied in P): 


in + l){xl 



Xn+l + Xn/e \ 
E -\- XfiXn-\-l ) 


n(l - xl/e^) 


Xu —1 T EX.jr 
£ “ 1 “ Xu — iXf 


= -yxn- 


- 1 


( 36 ) 


where e = e®“. This equation allows to represent Xn+i as a function of n, Xn-i and Xn'- Xn+i = 
$(n,x„_i,Xn). <h(n,u,u) maps the torus T'^ = x = {(u,u) £ C : |u| = |u| = 1 } into 
and has the following properties: 

• Vn £ N it is a continuous map on Aj x Aj where Aj = {e*^ : /3 £ [ 0 , a]}. 

• For (tt, u) £ A/ X Ai holds true <l>(n, u, v) € AjU Ajj U A/y, where A// = {e®^ : /3 £ (a, tt]} 
and Ajv = {e®^ :/?£[« — 7r,0)}, i.e. x can not jump in one step from Aj into Am = 
{e®^ : j3 £ (—TT, a — vr)}. 

These properties guarantees that there exists the unitary solution Xn = e®“" of this equation 
with xq = in the sector 0 < «„ < ct- This solution corresponds to embedded Z'^. Equation 

El is a special case of more general reduction of cross-ratio equation (see MM for the detail). 
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5 Circle patterns and Log 

For 7 = 2 formula m gives infinite Ri+i- The way around this difficulty is re-normalization 
/ —> (2 — 7)7/7 and limit procedure 7 —> 2 — 0, which leads to the re-normalization of initial 
data (see | 3 )- As follows from m this re-normalization gives: 

Ro = 0, Ri+i = 1, Ri = (37) 

a 


Definition 5 Circle pattern has radius function specified by Q),0) with initial data 
The symmetry 

7^2-7. 

of 0 )® stated in Lemma His the duality transformation (see (Hj). Smooth analog / 
holomorphic functions f{w),f*{w) is: 

df{w) df*{w) ^ ^ 
dw dw 

Note that log*(ui) = w‘^/ 2 . 

Definition 6 Circle pattern Log is a circle pattern dual to . 

Theorem 2 Discrete maps corresponding to circle patterns Z’^ and Log are embedded. 

Proof: The circle radii for and Log are subject to 0,® with 7 = 2 and 7 = 0 respectively. 
For these values of 7 Proposition|l]is true: the proof is the same since and Log are immersed. 
Due to Lemma 121 it suffices to prove the property m only for Z^. 

Consider the discrete conformal map ^^Z'^ with 0 < 7 < 2 . The corresponding solution Rz 
of 0,® is a continuous function of 7. So there is a limit as 7 ^ 2 — 0 , of this solution with 
the property m, which is violated only for z = 0 since i?o = 0- 
Q.E.D. 


(38) 
r for 


6 Concluding remarks 

Further generalizations of discrete Z'^ and Log are possible. 

I. One can relax the unitary condition for cross-ratios and consider solutions to 

Q{fn,m:fn+l,m}fn+l,m+l}fn,m+l') = K C (39) 

subjected to the same constraint 0 with the initial data 

„i-ya 

/i,o = 1, /o,i =-• (40) 

Kj 

This solution is a discrete analog of Z"^ defined on the vertices of regular parallelogram lattice 
(see Fig. |^. However, thus obtained mappings are deprived of geometrical flavor as they do not 
define circle patterns. 


14 





Figure 6: Discrete k = 2, a = 7r/2. 


II. Another possibility is to de-regularize prescribed combinatorics by projection of on a 
plane as follows (see [12] )■ Consider C R". For each coordinate vector Oi = 
where ej = Sj define unit vector in C = so that for any pair of indexes i,j vectors 
form a basis in R^. Let D £ R"^ be some 2-dimensional connected simply connected cell complex 
with vertices in Z”. Suppose 0 G D. (We denote by the same symbol D the set of the complex 
vertices.) Define the map P : D —> C by the following conditions: 

1) P(0) = 0, 

2) if X, y are vertices of D and y = x + ei then P{y) = P{x) + 

It is easy to see that P is correctly defined and unique. 

We call D a projectable cell complex iff its image u) = P(D) is embedded, i.e. intersections of 
images of different cells of D do not have inner parts. Using projectable cell complexes one can 
obtain not only regular square grid and hexagonal combinatorics but more complex ones, i.e. 
combinatorics of Penrose tilings. 

It is natural to define discrete conformal map on a; as a discrete complex immersion function 
/ on vertices of u preserving cross-ratios of w-cells. The argument of / can be labeled by the 
vertices x of D. Hence for any cell of D, constructed on e^, ej the function / satisfies the following 
equation for cross-ratios: 


Qifxi fx+e]^i fx+ei^+eji fx+Sj) — 6 




(41) 


where is the angle between and taken positive if (^kj^j) has positive orientation and 
taken negative otherwise. Now suppose that / is a solution to m defined on the whole Z^. 

Conjecture 1. Equation I2IP is compatible with the constraint 


ifx = 


S=1 


ifx+es - fx){fx - fx -es) 

/x+e. - /x-e. 


(42) 


For n = 3 this conjecture was proven in |1|. 

Now we can dehne discrete Z'^ : w ^ C for projectable D as solution to restricted 

on D. Initial conditions for this solution are of the form © so that the restrictions of / on 
each two-dimensional coordinate lattices is an immersion defining circle pattern with prescribed 
intersection angles. 
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This definition naturally generalize definition of discrete Z'^ given in [7j for Q = {(k,l,m) : 
k + I + m = 0, ±1}. 

Conjecture 2. Discrete : w —> C is an immersion. 

Schramm m showed that there is square grid circle patterns mimicking Erf(\/iz) but an analog 
of Erf (z) does not exist. The obstacle is purely combinatorial. There is a hope that combinatorics 
of projectable cells can give more examples of discrete analogs of classical functions. 
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